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The scattering of conduction electrons off neutral donors depends sensitively on the relative orien- 
tation of their spin states. We present a theory of spin-dependent scattering in the two dimensional 
electron gas (2DEG) of field effect transistors. Our theory shows that the scattering mechanism 
is dominated by virtual transitions to negatively ionized donor levels. This effect translates into 
a source-drain current that always gets reduced when donor spins are at resonance with a strong 
microwave field. We propose a model for donor impurities interacting with conduction electrons in a 
silicon transistor, and compare our explicit numerical calculations to electrically detected magnetic 
resonance (EDMR) experiments. Remarkably, we show that EDMR is optimal for donors placed 
into a sweet spot located at a narrow depth window quite far from the 2DEG interface. This allows 
significant optimization of spin signal intensity for the minimal number of donors placed into the 
sweet spot, enabling the development of single spin readout devices. Our theory reveals an inter- 
esting dependence on conduction electron spin polarization p c . As p c increases upon spin injection, 
the EDMR amplitude first increases as p 2 c , and then saturates when a polarization threshold pr is 
reached. These results show that it is possible to use EDMR as an in-situ probe of carrier spin 
polarization in silicon and other materials with weak spin-orbit coupling. 

PACS numbers: 72.25.Dc; 76.30.-v. 85.75.-d 
I. INTRODUCTION 

FM source ] gate | drain 

Electron spins carry much promise as state variables 
for scaled classical logioi and as quantum bits (qubits) 
in quantum computer architectures^ A key challenge of 
current research in spin-based electronics ( "spintronics" ) 
and quantum computation is to devise methods to probe 
spin polarization in semiconductors with weak spin-orbit 
coupling, 3 like silicon and silicon-germanium alloys. Sil- 
icon combines many special features that make its elec- 
tronic spin a promising basis for classical and quantum 
logic devices: Band structure properties such as weak 
spin-orbit coupling and indirect band gap, and the pos- 
sibility of preparing a nuclear spin free environment lead 
to extremely long spin relaxation times in comparison to 
other semiconductors." 

However, the same features that lead to long spin re- 
laxation times also make spin detection more challeng- 
ing. Optical methods for spin detection such as Faraday 
and Kerr rotation have been used successfully in III-V 
semiconductors^^ but unfortunately, these methods are 
inapplicable to silicon and related materials, whose spin- 
selective optical transitions are extremely weak and inef- 
fective. 

Recently, injection and detection of conduction elec- 
tron spin polarization was demonstrated in silicon. In 
Ref. Q, a Si device with a ferromagnetic (FM) source 
was interfaced with GaAs to enable optical spin detec- 
tion, while in [| a spin valve based on a FM-Si interface 
was used. 

Both experiments achieve spin detection at the in- 
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FIG. 1: (Color online) Silicon transistor with a ferromagnetic 
source, inspired by the devices in 001 • Our theory shows that 
doping with different types of group V donor species along the 
transistor channel enables the determination of the conduc- 
tion electron spin polarization distribution along the channel. 
Each donor species has a distinct set of hyperfine-split elec- 
tron spin resonance frequencies, whose amplitude is directly 
proportional to the local carrier spin and donor electron spin 
polarizations. 



terface, but no local or "in-situ" spin detection. The 
interface limits spatial resolution and introduces addi- 
tional scattering effects. In studies with accumulation 
channel silicon transistors, spin-dependent scattering was 
demonstrated by electrical detection of magnetic reso- 
nance (EDMR)^ii£ 

In this article we present a theory of spin-dependent 
transport in the scattering of spin polarized conduction 
electrons in the two dimensional electron gas of field effect 
transistors. Based on the scattering theory, we propose 
a donor-based approach to the problem of in-situ carrier 
spin detection in silicon. At the limit of a single donor 



impurity, this mechanism also provides a path for the 
readout of single donor spin states that does not require 
proximity to (or presence of) charge traps . 2 ] 11 i 12 i 14 

Consider a silicon transistor with a series of different 
group V donor impurities implanted along its conduc- 
tion channel, Fig. [TJ Spatially resolved characterization 
of carrier spin polarization is obtained by measuring the 
source-drain current when each donor species is in reso- 
nance with a microwave field. 

This is possible due to spin-dependent scattering 
(SDS), which relies solely on the symmetry of impu- 
rity scattering events3 i 20 ' 21 While this effect was believed 
to be observed through EDMR in Q amidst competing 
heating effects, it was only recently that SDS was con- 
firmed experimentally using a small ensemble (6 x 10 6 ) of 
implanted antimony (Sb) donorsJ^ In the following, we 
present a general microscopic theory of SDS, its efficiency 
and tunability in a two dimensional electron gas (2DEG) . 
We show that the core physical mechanisms underlying 
SDS are virtual transitions mediated by negatively ion- 
ized donor levels. Our results apply to all semiconduc- 
tors, but are particularly useful as a probe of spin po- 
larization in group IV and related materials with weak 
spin-orbit coupling, whose spin detection is intrinsically 
challenging. Our physical optimization of this effect has 
exciting consequences for spintronics and single donor 
electron or nuclear spin readout f22, 

We remark that SDS is quite different from other 
EDMR experiments based on spin-dependent recombi- 
nation of electron-hole pairs i 13 ' 14 ^ 15 ' 16 ! 17 ' 18 ^ 19 As showed 
by Kaplan, Solomon, and Mottj^ spin-dependent recom- 
bination has important contributions that do not depend 
on carrier spin polarization p c . Hence it is difficult to use 
recombination based EDMR as a probe of carrier spin 
polarization p c . 

This article is organized as follows. Section HT1 presents 
a general theory of SDS and its detection by EDMR in 
the case of thermal equilibrium in an external magnetic 
field (no spin injection). Section IIIII describes our ex- 
plicit numerical calculations, together with a comparison 
between theory and two experiments^^ and a discussion 
on how SDS can be optimized with respect to donor im- 
purity location. Section [TV] discusses the validity of our 
perturbation theory approach, by showing explicit calcu- 
lations of the Kondo temperature as a function of donor 
impurity location. Section [V] generalizes our theory to 
the case of non-equilibrium spin injection, and discusses 
how EDMR scales with carrier spin polarization. Sec- 
tion [VI] presents our conclusions. 



II. THEORY OF SPIN-DEPENDENT 
IMPURITY SCATTERING AND ELECTRICALLY 
DETECTED MAGNETIC RESONANCE 

Impurity scattering depends on whether the two- 
particle states formed by a conduction and an impurity 
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FIG. 2: (Color online) Band diagram along the z direction 
perpendicular to the 2DEG plane, as shown in Fig. [1] Here 
E c is the conduction band edge, ef is the 2DEG Fermi en- 
ergy, e a b,i is the ground state energy for the Z-th subband, 
and ei is the ground state energy of a neutral donor (singly 
occupied). We also show the donor impurity ground state 
wavefunction $1, and the conduction electron wavefunction 
9k,i with its characteristic length zq. As we show below, only 
donors located at R 3> zo can be detected by EDMR. 

electron are in a singlet (S) or triplet (T) configuration, 

*£/ T = ±= [* fc Kn)$i(r 2 ) ± * w (r a )$i(n)] . (l) 

Here ^>ki and $1 are orbital wavefunctions, respectively 
of a conduction electron with momentum k in the Z-th 
subband, and a localized donor impurity electron in the 
ground state. Fig. [2] shows the band diagram and poten- 
tial profile for conduction and impurity wavefunctions in 
the 2DEG. 

Access to additional channels for virtual scattering in 
the singlet state (Fig. [3]) translates into distinct neu- 
tral impurity scattering times for triplets vs. singlets, 
t t 7^ t~s- With spin polarizations p = (pj — pj.)/(pf +Pi), 
for conduction (p c ) and impurity (pi) electrons, the occu- 
pation fraction for singlets is given by ps — (1 — Pip c )/4, 
and the fraction for triplets pr — 1 — Ps = (3 + pip c )/4. 
Due to this difference in singlet vs. triplet scattering 
times, the device current is directly related to the spin 
polarizations according to 

I oc ( r ) = (p S Ts + Pttt) oc J (1 + ap c Pi) ■ (2) 

Here we introduced the parameter a as a figure of merit 
for the SDS effect. Eq. © forms the basis of the SDS 
mechanism of EDMR detection, which only occurs when 
the carrier spin polarization p c is non-zero. At ther- 
mal equilibrium, and in the low temperature "degenerate 
limit" (ksT <C Minijep — e s b,l}) 5 Pc arises due to Pauli 
paramagnetism, 

p c = "T-"i = WE . (3) 

n ] +n l 2 [e F - |(e sb ,i + e sb , 2 )] ' 
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FIG. 3: (Color online) Contributions to neutral impurity scat- 
tering, (a) The largest individual contribution to singlet scat- 
tering is a virtual transition to a negatively ionized donor 
level, where two electrons occupy the same (ground) orbital 
donor level e%. Such a transition is symmetry forbidden for 
triplets due to Pauli exclusion, (b) Triplet scattering occurs 
due to virtual transitions into negatively ionized donor levels 
where one electron occupies the level ei, and the other elec- 
tron occupies an excited orbital state ej > e\. These same 
transitions into excited orbital states are also allowed for sin- 
glets. 



Here n a , a =|, J, is the electron density for spin species 
cr, B is the external magnetic field, gp/Ti — 7 e = 
1.8 x 10 7 (sG) _1 is the free electron spin gyromagnetic 
ratio, and e s b,i, e s b.2 are the 2DEG subband ground state 
energies. 

Under microwave irradiation, the donor spin transition 
rates are given by W = W-\ — Wj. (the subscript refers 
to the initial state of the transition), 



W» = 27r|7eB±| 
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where B± is the amplitude of the microwave field, u) m — 
7e(i? + h m ) is the resonance frequency of the donor with 
nuclear spin state m = —I, —I + 1, . . . , +1 (h m is a hy- 
perfine shift), and ^ is the inhomogeneous linewidth. 
The microwave rate W competes with the donor spin- 
flip rate -jr — Tj + T| in order to determine a steady 
state impurity spin polarization, 

PiO 



1 + 22W' 



(5) 



where p^ = tanh ( ® k ^ T J is the equilibrium donor spin 
polarization. 

When donor atoms are at resonance (uj — ui m ), and 
satisfy the saturation condition (-f e B±) 2 > l/(TiT|), 
their spin polarization vanishes, pi ~ 0, effectively reduc- 
ing the device current by A// Jo ~ ap c pio- In the "clean 
limit", when no other scattering sources are competing 
with donor impurity scattering, the figure of merit a is 
given by 



J2i ( 3t t,; + t s ,i) 



(6) 



where we added a subscript I to account for different 
subbands. Note that the saturation condition depends 



critically on the microwave amplitude B±, donor spin- 
flip rate Jr- and resonance linewidth 

In transport measurements, SDS is always in competi- 
tion with other scattering sources, such as surface rough- 
ness and lattice defects. We describe these phenomeno- 
logically by adding a rate — : — ► — H — — . Usually, 

a J J o To Ts/T TQ r S / T 

transport is dominated by scattering sources other than 
neutral donors (— » ), leading to a reduced figure 

V T T S/T ' 

of merit «~^E/f^ L 

4 i—'l \T S ,l Vr,i t 

The scattering amplitudes to second order in the Born 
approximation are given by 
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and the scattering rates due to a single donor impurity 
are 
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(8) 



where 6 k k' is the angle between incoming (fe) and outgo- 
ing (k 1 ) conduction electron wavevectors, with |fc'| = k 
averaged over the Fermi distribution. In the degenerate 
limit we may set k equal to the Fermi wavevector fcp, but 
a more general averaging procedure will be considered be- 
low. Here p = ^ is the 2 DEG energy density, with A 
the 2DEG area. The energy of the conduction electrons 

is assumed to be e H = e sb ,i + -5^, with fc ± = (k x ,k y ), 
m* the effective mass, and e s t>,i the ground state energy of 
the Z-th subband. Eqs. ([7|) and {8} do not depend explic- 
itly on electron Zeeman energy because we assumed that 
the g-factor of the donor is the same as the g-factor of 
the electrons [while it is essential to consider different g- 
factors for EDMR detection, e.g. in Eq. the g-factor 
differences give only a negligible effect in the calculation 
of transport properties, provided the external field B is 
not too large]. 

The second order contribution in Eq. ([7]) is a sum over 
virtual (intermediate) states, and is generally much larger 
than the 1st order contribution. Fig. [3] illustrates the 
most important channels for neutral impurity scatter- 
ing. These depend on the two-particle Hamiltonian Ti. 
only through the donor single-particle orbital energies e,, 
donor "on-site" Coulomb repulsion U, and the overlap 
between conduction and impurity wavefunctions. The 
largest individual contribution is a virtual transition to 
(ex + U), which is a negatively ionized donor level with 
two electrons in the ground orbital state [Fig. [Ha)] . This 
transition is only allowed for singlets due to Pauli exclu- 
sion. Since singlet scattering also has the channels for 
excited orbital states available, it is always stronger than 
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triplet scattering (which has access only to the excited 
states) and in general a > 0. 

In these neutral donor scattering events, the total spin 
of conduction and donor electrons is always conserved, 
but their respective spin states may be exchangedi 23 i 24 
In our notation, T a is the rate for a donor with spin a 
to exchange spin with a conduction electron having spin 
—a; using a = +, — to denote the spin state, the rates 
can be written as 

2tt 



and C ee is the Coulomb electron-electron repulsion, 



r CT — -r - ^ \Jkl,k'l' | 2 S(tkl — Cfc'Z') 
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Here Jfcj.fc/^ - -i W>fc ,j, 
amplitudes, and 



Aui uni are exchange scattering 



f(e) = 



(10) 
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are Fermi functions. After converting the sums into in- 
tegrals over k space, the donor spin-flip scattering rate 
becomes 
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(11) 

As a check, note that the equilibrium donor spin polariza- 
tion is given by pio = (Fx— rj)/(Fj+Fj) = tanh 

When UbT 3> gpB, ^- scales linearly with the tempera- 
ture, while in the opposite regime it scales linearly with 
B. 

When other donor spin-flip mechanisms are active, we 
must add their rate to Eq. |9]). Nevertheless, we will see 
below that exchange scattering is the dominant contri- 
bution for donors in silicon in a wide parameter range. 



III. MODEL CALCULATIONS AND 
COMPARISON TO EXPERIMENT 

A. Model Hamiltonian and virtual two-particle 
donor states 

The two-particle Hamiltonian TL that models the cou- 
pling between a conduction plus a donor impurity elec- 
tron can be written explicitly as 

H = H (l)+H (2)+C R (1)+C R (2)+C ee (l, 2)-2e F , (12) 

where for notational simplicity (1) denotes coordinate T\, 
and (2) denotes coordinate r 2 . Here Ho = T+V z (z) is the 
translational invariant single-particle Hamiltonian, with 
T the kinetic energy and V z (z) the 2DEG confinement 
(triangular at low z but flattens out at high z, see Fig. [5]). 
C R is the attractive Coulomb potential of the impurity, 



C R (i) 
C ee (l,2) 



e 2 e -g T pk<-R| 

k |n - R\ 

g2 g-«TF|r-i-7- 2 | 

k \r 1 -r 2 \ 



(13a) 
(13b) 



where R = Rz is the location of the donor impurity, 
and k w 12 the dielectric constant for silicon. Since we 
are in the degenerate limit, it is important to account for 
screening; we use Thomas-Fermi screening with wavevec- 
tor 
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(14) 



where n 2 D is the 2DEG density and cf>(z) the subband 
wavefunction defined below. Note that the last term in 
Eq. (fT2)) is the chemical potential times the number of 
particles. 

In order to compute the scattering amplitudes 
[Eq. ([7])], we need to choose a set of virtual states \i s ' T ) 
that forms a complete basis for the two-particle Hilbcrt 
space. An important insight is that the conduction elec- 
tron may hop into the impurity site, form a negatively 
ionized virtual donor state, and then hop back into the 
Fermi sea. This motivates the choice of a molecular or- 
bital basis of negatively ionized donor states: 



iS/T) = 71 [l$1<l5 ' i) ± l$l$1>] ' z = 1 ' 2 ' 3 '- 



(15) 



where $i(r) is a donor orbital state with single particle 
energy e^, satisfying (Tt +Cr) = £i|$i). Note that 
i = l refers to the ground orbital state, and i = 2,3,... 
refers to other excited donor single-electron orbitals. 
From Eq. (fT5"|) we see that the state |1) only exists as 
a singlet, |1 T ) = 0. 

The simplest molecular orbital approximation is to as- 
sume 



n\i s ^ T )^(e 1 + e l + U-2e F )\ i s / T ). 



(16) 



Using this relation, we see that the 2nd order contribu- 
tion in Eq. ([7]) depends only on the energies e^, U, and 
on the overlap integral = (^kil&i)- Note also that 
states of the type \$i& 3 ) ± with i,j ' ^ 1 do not 

contribute to Eq. ([7J [they are orthogonal to Eq. (fT])] . 



B. Microscopic model for single particle states and 
overlap integral 

Bulk silicon has a six-fold degenerate conduction band, 
with energy minima located at k = 0.85 x 2vr/(5.43 A) 
along the set of (100) directions in the Brillouin zone. 
Its effective mass is anisotropic. Each valley has a heavy 
mass equal to 0.98m e along the valley direction, and a 
much lighter mass m* = 0.19m e in the perpendicular 
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plane (m e is the free electron mass). In the presence of 
a (001) interface, electrons in the valleys along ±(001) 
will have considerably lower energy a consequence of 
the heavy longitudinal massi 2 ^ This leads to the fol- 
lowing two-subband model for the conduction electron 
wavefunctions: ^kl(f) = ipi(z)e' lk " r± /y/A, with A the 
2DEG area and I — 1,2 a subband label. The subband 
wavefunctions are given by ip\ (z) = 0(z)\/2cos (koz) and 
ij)i{z) = (f){z)y2 sin (koz). For z > 0, we use the Takada- 
Uemura envelope function, 

*> = \/5(!;) e " K * r ' (17 > 

that is known to be an excellent analytic approximation 
to the self-consistent 2DEG ground state (see p. 468 
of Ref. [H). We assume <f>(z) = for z < 0. The 
characteristic length scale Zq models the 2DEG width, 
that can be controlled electrically by the gate voltage. 
Our choice of ipi corresponds to pure imaginary inter- 
valley scattering at the interface (other choices of inter- 
valley scattering phase simply give a phase shift to the 
cosine and sine functions, that is equivalent to changing 
the position of the interface) The conduction electron 
states are assumed to satisfy Tio\^k,i} = ^kl\^k,l), where 
€ki = 2nF^i ~^ e st>,i- Each subband has energy e s b.i at 
k = 0, with valley-splitting (e sbi 2 — e s b,i)/&B =0 — 100 K 
depending on interface quality and device geometry . 28 ! 29 
We now discuss model wavefunctions for the donor im- 
purity. We expect that donors located too close to the 
interface will have quite high ^- , and the EDMR satura- 
tion condition (•y e B±) 2 > l/(TiT 2 *) will not be satisfied. 
For this reason, we expect that only donors located far 
enough from the interface (at depths R 3> zq) can con- 
tribute to EDMR. Therefore it is sufficient to use donor 
wave functions that are good approximations in the bulk. 
Those are the Kohn-Luttinger wavefunctions ~ 

6 

$,(r - R) = F^r - R)J2 '/,,<'"' <r " • (18) 

For simplicity, we used an isotropic envelope (r — R) w 
e~ r '/< /y/iraf, with r' 2 =x 2 +y 2 + (z- R) 2 . The Bohr 
radius depends on the hydrogenic principal quantum 
number n according to a* w n 2 a* (we used a\ — 18.6 A). 
The principal quantum number n relates to the subscript 
i according to n = (i div 6) + 1. Below we will list the 
donor energy levels with the zero of energy chosen at flat 
band (the bulk conduction band edge). 

The valley vectors kj are given by fei = +k x, k 2 = 
—kox, k^ = +koy, k^ = —kgy, k$ = +koz, and 
kg = —koz. The corresponding rjij are conveniently writ- 
ten as 6-dimensional vectors. Each hydrogenic envelope 
has six different valley-split states, classified by symme- 
try. The lowest energy state for Sb has t\ = —45 meV 
with Al symmetry, tjai — 7^(1) 1) 1) 1j h l)^This is fol- 
lowed by a three-fold degenerate Tl symmetry level, with 



energy e 2 = —33 meV and tjti = ^(1,-1,0,0,0,0), 
^(0,0,1,-1,0,0), ^(0,0,0,0,1,-1). Note that only 
the last state couples to the (001) subbands. The highest 
energy level has E symmetry with energy £3 = —31 meV 
and is two-fold degenerate, tje = |(1, 1,-1,-1,0,0), or 
i(l, 1,0,0, —1, —1). Again, only one of these couples to 
electrons in the lowest energy subband. For principal 
quantum number n > 2 we assumed e„ = ei/n 2 with 
degenerecence 6n 2 . Finally, we used U = 43 meV for the 
on-site Coulomb repulsion. This corresponds to a bind- 
ing energy of Ed- = — (ei + U) = 2 meV; this is quite 
similar to the binding energy of 1.6 meV measured for 
phosphorous (P) in silicon, 32 and calculated by Oliveira 
and Falicov.-- We are not aware of measurements of Eo- 
for Sb impurities. 

One important point is that the overlap integral be- 
tween conduction electrons in the l-th. valley and the im- 
purity electron Su,i = (^u\^i) is strongly oscillatory 
on R. For notational convenience, define s% — cos(/co-R), 
and S2 = sin(fco-R), and set S3 = s\. After a simple 
calculation we see that Ski,i oc s; apart from a smooth 
envelope whenever the valley symmetry of donor orbital 
i is "even" (oc cos [ko(z — R)]), and Ski,i oc s/+i when the 
valley symmetry of orbital i is "odd" (oc sin [ko(z — R)])- 
These oscillations occur due to valley interference, in a 
similar fashion as the exchange oscillations between two 
donor impurities in silicon. 34 ! 35 ! 36 

C. Contributions to the scattering amplitude 

We start by evaluating explicitly the first order contri- 
butions to Eq. ([7]) . To get an idea of the order of magni- 
tude of each contribution, we will quote numerical values 
for a single donor located at R — 6zq, with parameters 
A = (lfim) 2 , e F = -9.32 meV, e sb ,i = e sb ,2 = -10 meV, 
and zq = 40 A. In this case the overlap integral between 
conduction electron at ep and donor ground state was 
Sfe.isi, with S kA = 1.31 x 10- 4 . 

Let's consider the first order contribution to Eq. 0. 
Using Wo|*«) = ew|*«) and [H +C fl ]|$) = ei|$), we 
get 

(*f/ T |W|*f/J) = (*w|CR|*fc/i/) + <*w$|C ee |* fc / r $) 
± [(ei - e F )Sk,iS^, A sisi' 

+($*ki|Ci2|**'j'*)]. (19) 

The first contribution on the right hand side is attractive 
Coulomb potential scattering (always negative), and the 
second contribution is repulsive scattering from the donor 
electron cloud. This has the same order of magnitude as 
the first term, but the opposite sign. Assuming s/ = 
si> = 1 (maximum value) these terms are respectively 
—8.1 neV, and +8.5 neV. The first term inside the brack- 
ets equals —0.61 neV, the second —0.3 neV, and the third 
+0.2 neV. Hence we have (^fi\H\^f, lt ) = -0.3 neV, and 
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O^jwl^l^fer) = +1-1 neV. Interestingly, in this partic- 
ular case the first order contribution to Eq. Ffy favors 
triplet scattering. Actually, computations for different 
donor positions R shows that the first order contribu- 
tion changes sign several times as the donor position R 
is varied. 

We now turn to second order contributions to Eq. 1(7} . 
We remark that these are always negative for singlets 
and for triplets because ei < e F < e\ + U < ei + U 
for i = 2, 3, . . .. The largest individual contribution is 
a virtual transition to a negatively ionized donor state 
with two electrons occupying the donor ground state. 
This contributes exclusively to singlet scattering. Using 
Eq. (fT6|) we get 



-2(2^ + U - 2e F ) 2 \S ktl \ 2 aw 
ei + U — e F 



(* s k i\n\i s )(i s \H\* s kll ,) 

eu —ei — U 

which equals —3.7 neV for the particular case consid- 
ered. This is more than 3 times larger than the 1st order 
contribution mentioned above. We checked several other 
parameter regimes, and found Eq. (|2"D|) to be 2-10 times 
larger than the 1st order contribution (singlet or triplet). 

All the other second order terms contribute equally to 
singlet and triplet. The largest of these involve virtual 
excited states to the n = 1 state with Tl or E symmetry. 
The Tl state contribution is 



_ g (ei + £4 + U - 2e F ) 2 \S k ,\\ Si +1 sy + i 
e^ + U - e F 

and the E contribution is 

3 (ei + e 6 + U - 2e F ) 2 |ff M | 2 sjsy _ 
2 e 6 + U - e F 



■ = -0.95 neV, 
(21) 

-0.47 neV. (22) 



The contributions with principal quantum number n > 
2 may be lumped together in a sum: 



^ 2 (ei +e n + U -2e F ) 2 , „ |2 



n=2 



U-e F 



\Sk,n\ 



x I 2 S * S *' + 3s i+ isp+i 
-1.3 x 10 4 neV, 



(23) 



since there are 3n 2 orbitals coupling to the 2DEG valleys, 
each with energy e n — ei/n 2 . The envelope of the overlap 
integrals is defined as 



Sk.i 



1 

71 



d z r(j>{z)- 



Akr± 



■F n {r-R). (24) 



We get A 2 = 390 neV, A 3 = 6170 neV, A 4 = 4164 neV, 
A 5 = 1367 neV, A 6 = 425 neV, A 8 = 50 neV. The A n 
reach a maximum for n = 3, and then decrease appre- 
ciably with increasing n because the Sk n 's become expo- 
nentially small. As a result, we get a good approximation 
by evaluating the sum up to n roax = 4. We found that 



the contributions for n > 5 produce negligible changes to 
our final result. 

Hence the sum of the second order contributions in 
Eq. ||7J) can be as much as 10 4 times larger than the 
first order contributions. It is interesting to note that 
the difference between singlet and triplet rates is de- 
termined by the difference between squared amplitudes: 
4 ~ ^: « (A% wv ) 2 - (Al hkn ,) 2 , which is of the order of 
10 4 neV 2 . This is several orders of magnitude larger than 
the square of the exchange scattering amplitude, that de- 
termines exchange scattering: jr oc (A^ k >i'~Afi k >i') 2 is 
only about ~ 10 neV 2 . Therefore, we see that exchange 
scattering is quite different than SDS, in the sense that 
the latter benefits from a large number of additional elec- 
tronic transitions. 



D. Dependence of EDMR parameters on donor 
depth 

We now show explicit numerical calculations of the pa- 
rameters determining EDMR detection, and discuss their 
dependence with donor depth R. The figure of merit a is 
shown in Fig. [4] and the donor spin flip rate y- due to the 
exchange scattering mechanism is shown in Fig. [5j Both 
a and y decrease appreciably as R increases relative to 
the 2DEG thickness z , and are quite sensitive to e F — e s b,i 
(or, equivalently, the 2DEG area density). We used 
e F — e s b,i = 0.68 meV, and considered two different cases: 
Valley degenerate with e s b,i = e s b,2 = —25 meV, and 
valley-split with e s b,2 = — 15 meV and e s b,i = —25 meV. 
Both a and are independent of the 2DEG area A; a 
depends sensitively on the scattering time due to other 
sources; we used tq = 0.4 ps which is equivalent to a 
typical transistor mobility of 4 x 10 3 cm 2 /(Vs) (^r is in- 
dependent of ro). We assumed T — 5 K and B — 0.36 T 
(a does not depend on B, and is nearly independent of 
T in the valley degenerate case). 

Fig. [5] shows that exchange scattering is several orders 
of magnitude stronger than conventional spin-phonon 
coupling of isolated donor impurities. For example, Fe- 
her and Gere 3 ^ measured -^r ~ 10 _1 s _1 for phosphorous 
donors in bulk silicon at T — 5 K and B — 0.36 T. 

Interestingly, the EDMR amplitude depends on 
whether the device temperature is larger than the valley- 
splitting energy or not. At large temperatures (or small 
valley splittings), a and y- are smooth functions of R 
(apart from tiny oscillations), because contributions from 
the two subbands complement each other (dark curves). 
For the opposite regime of temperature lower than valley- 
splitting, the parameters are strongly oscillatory on R 
(grey curves). In this case tjt becomes quite small for 
some donor positions close to the interface, suggesting 
that a fraction of the donors located at "lucky sites" 
might be detectable by EDMR. In Figs. 0] and [5] we used 
filled circles to denote the actual silicon lattice sites that 
may be occupied by substitutional donors. 
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FIG. 4: Figure of merit a as a function of donor depth 
R divided by the 2DEG thickness Zq. Dark curve: Val- 
ley degenerate case. Grey curve: Valley split case, with 
e s b,2 — £sb,l S> ksT. Filled circles denote actual sites of the 
diamond lattice, which may be occupied by donors. 
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FIG. 5: Donor impurity exchange scattering rate as a 
function of donor depth R/zo (same parameters as in Fig. [4j . 
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TABLE I: Comparison between theory and experiment at 
T — 5 K and B = 0.36 T with an inhomogeneously broad- 
ened linewidth of — = 2 G. Our theoretical calculations 

of A7/7o were obtained by multiplying the donor density by 
the scattering rate due to a single donor, and integrating over 
donor depth. We used 2DEG thickness zo = 50 A, valley- 
degenerate subband energy e s b,i = e s b,2 = —25 meV (zero en- 
ergy is at flat band), and transistor mobility 4 x 10 3 cm 2 /Vs 
(corresponds to tq = 0.4 ps). 



E. Comparison to experiment 

Our theory can be directly compared to experiments. 
Ghosh and Silsbee- measured AI/I = (I a s res. — 
^at ros.)/^o ~ — 10 -5 at high power (the current increased 
upon resonance), while at lower power they found that 
AI changed sign. In this study, a silicon transistor bulk 
doped with phosphorous impurities was used. Lo et al^ 
measured much lower amplitudes | AI\/Iq = 10~ 8 — 1CP 7 , 
with antimony donors implanted only into the transistor 
channel. Unfortunately the use of derivative detection 
did not allow the determination of the sign of AI. 

Our theory disagrees in sign with the measurements of 
@, suggesting that their high power signal was not due 
to the SDS mechanism, instead it was likely due to the 
2DEG heating mechanism. 38 

We now compare our explicit numerical results with 
the experimental data of Lo et al. 10 Our theory al- 
lows the prediction of the device current for a given 
donor distribution without any fitting parameters. De- 
note Nj = AAz rid(zj) the number of donor impuri- 
ties located in the interval Az = — Zj, with n^{zj) 
the volume density of donors in the Az-thick layer. The 
EDMR signal averaged over the donor profile is simply 
given by 

— p-^ =PcY\Nja{zj)pi{zj). (25) 
Jo 

3 

We obtained the depth dependent donor density nd(z) 
from Secondary Ion Mass Spectroscopy measurements 
(SIMS). Our explicit model calculations of the EDMR 
current are shown in Table [J for all measurements made 
in The EDMR response predicted by theory is in 

reasonable agreement with experiment. 

F. Physical optimization of EDMR 

Table U shows that the SDS contribution to EDMR 
is much weaker than anticipated on the basis of previous 
measurements.— We now describe a physical optimization 
that aims at maximizing | AI\/Iq for the minimal number 
of donor impurities. 

The SDS parameters a and ^- decrease rapidly with 
increasing donor depth R, an effect that has important 
implications for the optimization of the EDMR ampli- 
tude. Notably, for given B±, the EDMR signal will be 
maximal for donors satisfying ^-(R*) w T^^eB^) 2 . The 
location R* is the closest one to the interface that satis- 
fies the saturation condition (j e B±) 2 > 1/(T-lT%), 

Fig. [6] shows the EDMR amplitude as a function of 
donor location, for a single donor implanted in the center 
of a transistor of area 0.1(/^m) 2 . We assumed microwave 
amplitude B± = 0.3 G, 2DEG density 1.1 x 10 11 cm 2 , 
with other parameters as in Table I. 

Interestingly, the EDMR amplitude depends on 
whether the device temperature is larger than the valley- 
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Sweet spot R' 

i 




Donor depth R/z 

FIG. 6: EDMR amplitude at T = 5 K as a function of 
donor placement, for zero valley-splitting (dark curve) and 
for valley-splitting equal to 100 K (oscillating gray curve, with 
filled circles denoting actual donor sites of the diamond lat- 
tice) . 

splitting energy or not. At large temperatures, |A7|/7o 
is a smooth function of 77, because contributions from 
the two subbands complement each other. In this case 
a single donor placed at the sweet spot 77* w 5.5 zq is 
optimally detected by EDMR (dark curve in Fig. [6] with 
(sbi = £sb2 = —25 meV). For the opposite regime of tem- 
perature lower than valley-splitting, |A7|/7o is strongly 
oscillatory on 77, and donors closer to the interface can 
be detected as well (oscillating gray curve in Fig. [H with 
e s bi = —25 meV and e s b2 = — 15 meV). 

Fig.[S]shows that EDMR is able to detect a single donor 
spin implanted in a transistor of area 0.1(/xm) 2 . For typ- 
ical 7 ~ 5 juA, current modulations A7 ~ 0.1 pA are 
detectable with standard techniques, provided the donor 
is placed at the sweet spot 77* . A single donor resonance 
may be identified by the presence of only one hyperfine 
satellite line out of the full satellite spectrum with 21+1 
lines (for donor nuclear spin I), for measurements within 
the nuclear spin-flip time^ 

The EDMR amplitude is directly proportional to the 
donor area density per monolayer. Therefore, if 6 x 10 6 
Sb donors were to be placed exactly at the sweet spot 77* 
in a large transistor of area ~ 10 3 (/im) 2 , the signal would 
be |A7|/7q ~ 10~ 5 , two orders of magnitude higher than 

in gg. 




Donor depth R/zo 

FIG. 7: Kondo temperature as a function of donor depth 
R/za (same parameters as in Fig. [4). Note that the Kondo 
temperature is exponentially small for R > 3zo; this justifies 
our perturbative approach. 



the lowest subband energy. Eq. ([55} applies when e\ + 
U > £f and £i < e s b < £f- The first term in the square 
root is a characteristic particle excitation bandwidth for 
electrons tunneling into the donor, while the second term 
is a hole excitation bandwidth. 

Fig. [7] shows our calculated Kondo temperature as a 
function of donor depth. Here we see that for 77 > 3zo, 
we have T^ndo <C 0.1 K. Hence, even at the lowest tem- 
peratures achievable experimentally, there is no Kondo 
effect for donors located at 77 > 3zo (7 1 ^> Tkondo). 

In the valley degenerate regime, only donors around 
77 ~ 5.5zo can be detected by EDMR, as shown in Fig. [6] 
(donors closer to the interface can be detected at higher 
power). Hence we see that donors located in this region 
satisfy T 3> Txondo for the lowest temperatures achiev- 
able in the laboratory, and our perturbation theory ap- 
proach is justified. In the valley degenerate case, both 
the EDMR current as well as the Kondo temperature are 
strongly oscillatory with donor depth, and again it can 
be seen that all donors that are detectable by EDMR 
have extremely low TKondo- Note that the behavior of 
TKondo as a function of donor depth is qualitatively sim- 
ilar to exchange scattering; as a consequence, all donors 
with y~ low enough to be detectable by EDMR have a 
corresponding quite small Tkondo- 



IV. KONDO TEMPERATURE AND VALIDITY 
OF PERTURBATION THEORY 

Our approach is based on perturbation theory [Eq. 0]. 
This is known to be a good approximation only when 
the temperature is larger than the characteristic Kondo 
temperature TKondo- For our problem, the Kondo tem- 
perature may be written as2£ 

fcsT Ko ndo = y/{e 1 + U - e F ) (e F - e sh ) | Jp|e~pW, (26) 

where J = (Jki,k>i>) is an average exchange energy at the 
Fermi level (averaged over subbands), and e s b is taken as 



V. SPIN INJECTION 

We now discuss the spin injection regime. Under spin 
injection, each spin species have different quasi-Fermi en- 
ergies, e_F| ^ efi^ For k B T -C Min{e F o- — £sb,i}, the 
carrier polarization is well approximated by 

p = n T - n l _ 2 (e FT - e Fl + gfiB) 
«T+ n l J2i ( e F T + £fi - 2e sb ,i) ' 

Here n a = J pdef a (e) is the 2DEG area density for each 
spin subband. 

The impurity scattering times Tg/x must now be cal- 
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culated for each spin subband: 



( T S/Td) c 



J pde (e - e sb ,i + <J^f-^j 



dfo 



~S/T,kl 



(28) 

where we added a spin subscript to the scattering times 
and the Fermi functions. The energy integrals are from 
(esb,i ^ er 2 ^-) to oo. The current is now calculated as in 
Eq. @, using ps = (1 — Pi)/4 for the spin up subband, 
and ps = (1 +pi)/4 for the spin down subband, with 
Pt = 1 — ps in each case. In addition, each {ts/t)u 
must be multiplied by its corresponding density n a . The 
source-drain current becomes 

I-Iq = Ej n r (( T ^)t ~ ( r ^)t) ~ n l (( T Ti)j - (rsi}j) 

J o E;"T ( 3 ( t t;)t + ( t s;)t) + n l ( 3 ( r ^)| + (tsi)i) 

(29) 

Interestingly, under spin injection the scattering times 
are intertwined with the densities n a , and a figure of 
merit a independent of p c can not be defined. Note 
that when (ts/t)| = ( t s/t)i, Eq. ([29]) becomes AI/I n = 
ap c pi, as obtained previously. 

The donor spin transition rates are also modified. Fol- 
lowing from Eq. ©, T a oc J de/_ CT (e)[l - / CT (e - cx^)], 
leading to 

1 27T 



x coth 



2fc B T 



The steady state impurity polarization becomes 
R0 = tanhf ^-^+ g ^ ). 



2fc s T 



(30) 



(31) 



Fig. [8] shows EDMR lineshapes under spin injection, 
using T — 5 K, _Bj^ = 0.3 G, and doping profile and 
other parameters similar to |l0| . The 2DEG density is 
fixed at 2 x 10 11 cm -2 , and the conduction electron spin 
polarization is varied between 2% and 75%. At low p c , 
the EDMR amplitude scales as p c Pio oc p 2 ., since a re- 
mains unchanged and pio scales proportional to p c due 
to exchange scattering [Eq. (|31[)] . In this regime, EDMR 
can be used as a local probe of carrier spin polarization. 

However, as p c is increased beyond the threshold 



Pt 



Max{2fc B T, g^B} 

(ei?| + EFJ.) — (fsbl + £sb2) ' 



(32) 



the spin- flip rate starts increasing as ^- (x p c , because 
any donor spin antiparallel to the 2DE(j spins will relax 
rapidly [note that Eq. ([30]) becomes proportional to p c 
when p c > pt]- By virtue of these larger ^j-'s, EDMR 
will excite donors further away from the interface, tend- 
ing to decrease AI/Iq. Remarkably, this effect competes 
against the p\ scaling, saturating AI/I for p c » pr- 
While this limits EDMR as a probe of carrier spin for 
Pc Pt, it also demonstrates that EDMR detection is 
optimal at p c ~ pt- 




3.58 3.59 3.60 3.61 

B (10 3 G) 



3.62 



FIG. 8: EDMR lineshape under spin injection, for carrier spin 
polarization p c = 2, 6, 10, 20, 30%. Inset: EDMR amplitude 
at resonance, as a function of p c . 



VI. CONCLUSIONS 

In conclusion, we presented a microscopic theory of 
spin-dependent scattering in the interaction of conduc- 
tion electrons with neutral donor atoms. Our results 
are based on an effective mass approximation. More 
sophisticated approaches are likely to reduce the valley 
oscillations, 36 with no modification to our conclusions. 
The considered mechanism requires temperatures lower 
than the impurity binding energy (T < 100 K), but even 
higher temperatures may be achieved using deep level 
magnetic atoms or clusters. 

We showed that SDS is determined by virtual transi- 
tions into doubly occupied donor states. As a result, SDS 
always leads to a reduced current upon EDMR satura- 
tion, since a > 0. 

A recent experiment^ demonstrates that the EDMR 
current due to SDS is indeed reduced upon donor spin 
saturation (Fig. 2(c) in [39]), in agreement with our pro- 
posed virtual transition mechanism. 

Spin-dependent scattering detection is challenging due 
to competing heating effects^ which are directly propor- 
tional to the number of donors present. Our finding that 
SDS arises solely from impurity spins located within a 
narrow depth window with respect to the 2DEG shows 
a path to significant optimization of spin signal intensity 
for a minimal number of donors placed into this depth 
window, and underpins the development of single spin 
readout devices. 

Our theory shows how the EDMR amplitude will 
scale with carrier spin polarization in the regime of spin 
injectioni^i^ Therefore, the monitoring of donor electron 
spin resonances can be utilized for the spatially resolved 
characterization of conduction electron spin polarization, 
providing a sensitive probe for optimization of spin injec- 
tion and spin transport in semiconductors with indirect 
band gap and weak spin-orbit coupling. 
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